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CONTINUOUS PARAMETRIC MODEL FOR CIRCUIT SIMULATION 

TECHNICAL FIELD OF THE INVENTION 

This invention relates generally to the field of computer aided design and analysis, and 
more particularly to an improved system and method for computer-based circuit modeling. 
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BACKGROUND OF THE INVENTION 

In the synthesis and analysis of electronic circuits and systems, computer design tools are 
pervasive. Circuit analysis software, such as various embodiments of the public domain 
simulation engine SPICE, allow engineers to efficiently test complex circuit designs using 
computer-based models of the circuits and circuit elements, prior to building physical prototypes. 
As these computer-based tools evolve, all but final design verification will ultimately be 
performed using computer-based modeling. 

However, as anyone with experience in using computer aided design tools is well aware, 
these tools are still somewhat limited. For example, a significant problem with many SPICE 
models is that they are designed in a piecewise fashion with different modeling equations 
characterizing different regions of the model. For example, a transistor model may include a 
parametric model for a voltage transfer function that has a first mathematical representation for a 
sub-threshold region, a second mathematical representation for a linear region of operation and 
yet a third mathematical representation to describe the saturation region of the circuit. This 
modeling approach can present problems in the continuity of device parameters (such as 
capacitances) as well as in the continuity of the derivatives of the parametric models. 

Numerical techniques have been proposed to deal with the discontinuities which result 
from piecewise linear modeling. A common solution for such piecewise linear models is to 
create a smoothing function. The purpose of the smoothing function is to maintain the 
continuity of the circuit model over an expected range of values. 

A popular smoothing function, representing the effective voltage of an electrical 
component, is set forth below in Equation 1. 



V ejr =V 0 -±[(V 0 -V-S) + ^(V 0 -V-5) 2 +4.S.V 0 ] 
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Equation 1 exhibits the properties that 

(2a) 

V 0 ifV»V 0 
V ifV«V 0 

and 

V eff = 0 if V=0 and V 0 > 0 (2b) 

Property 2(b) is an important attribute of a smoothing function. For example, the 
smoothing function of Equation 1 can be used to model the drain-to-source voltage (V^) of a 
FET transistor, in the following manner: 

V dseff = V dsal - \ {(V^ - V ds -S) + JiV^-V^-Sy+A.S.V^ } (3) 

The smoothing function of equation 3 has the property that at = 0, the effective drain- 
to-source voltage (V dseff ) is also zero. Additionally, in accordance with the property given in Eq. 
(2a), approaches if is much greater that V^, but becomes itself if is smaller 
that V^. 

The constant 8 is used to provide a smooth transition at the corner point of V=V 0 . If V eff 
were defined as Eq.(l) except with 5=0, i.e., 

V eff = Vo \ {(To ~ V) + V(^0-n 2 } (4) 

then the resulting function exhibits a piecewise continuous nature and its first derivative would 
be discontinuous at the point V=V 0 . 

A typical value used for 5 is 0.02. For a certain case when V 0 =l , V efr obtained from 
Eq.(l) as a function of V is shown in Fig. 1. Figure 1 is a graph which illustrates that V efT is 
smooth and achieves both properties specified in Eq. (2). However, there are two severe 
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drawbacks to the smoothing function of Eq. (1). First, when V G assumes a negative value a 
discontinuity is encountered. Figure 2 is a graph which illustrates the same calculations used to 
generate the graph of Figure 1, except that V 0 is equal to -1 in this example. Figure 2 
demonstrates that even though V eff converges properly to V 0 when V + V 0 , a discontinuity 200 
5 occurs near the turning point. 

The second problem of Eq. (1) results when V 0 is identically zero. In this case, V eff 
calculated from Eq. (1) is shown in the graph of Figure 3. Figure 3 illustrates that V efr is 
piecewise continuous and, therefore, that the first derivative is discontinuous at V o =0. This 
problem results from the fact that if V o =0, the term 4$8$ V 0 inside the square root ceases 
10 providing the smoothing action because it is identically zero. 

^ When the models used by a simulation engine exhibit discontinuities, the resulting 

O 

H simulation may yield incorrect results, fail to converge or even encounter run-time computer 

\.n 

li errors. While known smoothing functions reduce the likelihood of such occurrences, care must 

p still be taken in determining the range of input values which are applied to many circuit element 

i ! * 15 models. Thus, there remains a need for an improved method of modeling component parameters 

n 

Q such that the aforementioned discontinuities are avoided and a broader range of input variables 

p can be applied in the resulting circuit models. 

P 

3 . 
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SUMMARY OF THE INVENTION 

Accordingly, a need has arisen in the art for an improved circuit models for use with 
simulation and analysis systems. 

In accordance with the present invention, a method of generating a continuous parametric 
model of an electronic circuit parameter having a base model includes the step of determining 
whether the base model exhibits at least one discontinuity over an allowable range of parameters. 
If the base model exhibits at least one discontinuity, 

at least one compensation function is applied to prevent the base model from exhibiting 
such discontinuities over the allowable range of parameters. The method also includes the step 
of determining whether the first derivative of the base model exhibits at least one discontinuity 
over the allowable range of parameters. If the first derivative of the base model exhibits at least 
one discontinuity, at least one compensation constant is applied to prevent a first derivative of 
the base model from exhibiting discontinuities over the permissible parametric range. 

Also in accordance with the present invention, a continuous parametric model of a 
physical circuit element is generated from a base model that defines a piecewise linear 
representation of the circuit element and exhibits at least one discontinuity over an allowable 
range of model parameters and/or exhibits a discontinuity in the first derivative of the allowable 
range of model parameters. At least one compensation function is included to prevent the 
continuous parametric model from exhibiting discontinuities over a permissible parametric range 
of positive and negative values. At least one compensation constant is included to prevent a first 
derivative of the continuous parametric model from exhibiting discontinuities over the 
permissible parametric range. 

Technical advantages of the present invention include providing improved component 
models for circuit simulators. In particular, the resulting continuous parametric models and their 
derivatives are free of discontinuities. 
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Other technical advantages of the present invention will be readily apparent to one skilled 
in the art from the following figures, descriptions, and claims. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present invention and advantages thereof, 
reference is now made to the following description taken in conjunction with the accompanying 
drawings, wherein like reference numerals represent like parts, in which: 

Figure 1 is a graph of a smoothing function for an exemplary circuit parameter (V eff ), as 
modeled in the prior art. 

Figure 2 is a graph of a smoothing function for an exemplary circuit parameter (V efT ), as 
modeled in the prior art. 

Figure 3 is a graph of a smoothing function for an exemplary circuit parameter (V efT ), as 
modeled using a smoothing function known in the prior art. 

Figure 4 is a block diagram of a circuit modeling and simulation system. 

Figure 5 is a flow chart illustrating a method of generating an enhanced continuous 
parametric model for a circuit element in accordance with the present invention. 

Figure 6 is a graph of an exemplary circuit parameter (V eff ) as modeled in accordance 
with the present invention. 

Figure 7 is a graph of an exemplary circuit parameter (V eff ) modeled with an enhanced 
continuous parametric model in accordance with the present invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

The preferred embodiments of the present invention and the attendant advantages thereof 
are best understood by referring now in more detail to the drawings in which like numerals refer 
to like parts. 

Figure 4 is a simplified block diagram illustrating an exemplary computer system which 
can be used to generate and use enhanced component models to model and simulate electronic 
circuits and systems in accordance with the present invention. The present invention is not 
specific to any particular computing platform, therefore, the computer system can take the form 
of a personal computer, workstation, network terminal or other computer devices. Generally, 
the system includes a processor section 410 which is operatively coupled to one or more 
computer readable mass storage devices such as disk drives, CD-rom drives, tape drives and the 
like. The processor 410 section can access the storage devices and import digital data 
representing simulation/analysis engine software 420 and one or more component model 
libraries 430. While shown in Figure 4 as separate structures, the program data for the 
simulation engine 420 and the data for the model library 430 can reside on a common storage 
device. The system will generally include an input device 440, preferably including a keyboard 
and digital pointer device such as a mouse, trackball, light pen, touch screen input device and the 
like. The system will also include a display device 450, such as a CRT screen. The input 
device 440, processor section 410 and display device 450 preferably cooperate to establish a 
graphical user interface (GUI). Of course, numerous variations on these basic components are 
possible without deviating from the present invention. 

In conventional circuit modeling, such as Compact7 and SPICE modeling, circuit 
components are described by models. The models, in turn, are generally a collection of 
mathematical representations of various device parameters which characterize the component, 
such as input/output transfer functions. Such mathematical representations are referred to herein 
as parametric models. Thus a particular component model is formed with one or more 
parametric models. In this disclosure, the term base parametric model refers to parametric 
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models, such as that described in connection with Eq. 1, which exhibits at least one discontinuity 
or the derivative of which exhibits at least one discontinuity. 

Referring to Figure 5, the method of generating an enhanced continuous parametric 
model begins with the acquisition of a base parametric model (step 505). The base parametric 
model is analyzed to determine whether the model is continuous over the desired full range of 
variables (step 510). If the base parametric model is not continuous, at least one compensation 
function is applied to eliminate the discontinuity (step 515). If the base parametric model is 
continuous, the method advances through logical connector 520 to determine whether the 
derivatives of the base parametric base model are continuous (step 525). If the derivatives are 
not continuous, a compensation constant is applied to the parametric model to eliminate the 
derivative discontinuity (step 530). After the derivatives have been found to be continuous, or 
have been compensated to be continuous, the method concludes with the result being an 
enhanced continuous parametric model (step 540). The enhanced continuous parametric model 
can be stored in model library 430 for use by the stimulation/analysis engine software 420. 

In accordance with steps outlined in Figure 5, a base parametric model is transformed 
into an enhanced continuous parametric model. As an example of the present invention, the base 
parametric model of Eq. 1 can be transformed into an enhanced continuous parametric model in 
accordance with the present invention. 

As noted in connection with Figure 2, Eq. 1 can exhibit a discontinuity when V 0 is less 
than zero. In removing the discontinuity in Equation 1, the term 8, a constant, can be replaced 
with a variable, 0, such as a mathematical equation. The term 0 can take the form of a number of 
functions which exhibit a small value when compared to V 0 . Equation 5, which is set forth 
below, illustrates this initial substitution. 




(5) 
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In Eq. (5), the function designated as 0 is referred to as a compensation function. A 
simple example of a compensation function which removes the initial discontinuities 
encountered in Eq. 1 is defined in the following manner: 

(6) 

f S ifV o >0 
eff \-S ifV o <0, 

where 5 is any arbitrarily small number, such as 0.02, as used in connection with Eq. 1. With 

this compensation function, the discontinuity problem illustrated in Fig. 2 is removed. However, 

the function of Eq. (6) is itself not continuous at V 0 = 9. An alternative compensation function 

can be defined as: 



io w)=y ( 7 ) 



where K is any positive number substantially greater than unity. Figure 6 is a graph of V eff of 
* Eq. (5) with 0 from Eq. (7) and K=100. As the graph of Fig. 6 illustrates, such an enhanced 

parametric model is satisfactory even when V 0 is negative. 

ST 

Z The enhanced parametric model set forth above in Equation 5, though representing an 

1 

« 15 improvement over the smoothing function of Eq. 1, still exhibits a discontinuity when V 0 is equal 

j 

I to zero. To correct for this effect, an enhanced parametric model of the present invention can be 

further defined as follows: 

v e ff = v *-\ fa -v-e) + J(y;-v-e) 2 + 4er;\ (8) 

where V 0 * includes a second compensation function to modify the terms with a V 0 component as 
20 follows: 

r;=F 0 +A.exp(-r 0 2 ) (9) 
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In this manner, V 0 * does not exhibit a value equal to zero even when V 0 is equal to zero. 
The constant parameters are preferably selected to be small, such that V 0 * approaches V 0 when 
V 0 is large, as set forth by the requirement of Eq. 2b. It will be appreciated that the expression 
for V 0 * given Eq. (8) is merely exemplary and that any number of specific functions exhibiting 
5 the aforementioned properties will suffice, thus, satisfying step 515. 

While the enhanced parametric model of Eq. (8) is substantially continuous, this function 
is only piecewise continuous when V o =0. Therefore, the derivative of this function will not be 
continuous at this point (step 525). In this case, the cause is that at V o =0, 4$ V o $0 is identically 
zero. To avoid this situation, a compensation constant, A, can be applied to the square root term 
10 of Eq. (8) in addition to the terms 4$V O $0 and (V 0 -V-6) 2 . 

1:3 

^ ^ Applying a mathematical identity, 

'■2 (F o -0) 2 +4#0#F o +2«K o #A + A 2 = (F o +0 + A) 2 , (10) 

:=3 an enhanced parametric base model can be defined as: 

I!-- V eff =V o -^{(V o -V-0-A) + ^V o -V-0) 2 +40V o +2V o A + 20*A + A 2 \ (11) 

i c 

t: 

O 15 With *=0.01, for example, the terms inside the square root other than (V o -V-0) 2 will never be 

□ zero even though V o =0. Recognizing that the overall value within the square root can be 

r ~ negative when the value of Vo is very small, and that the property dictated by Eq. (2a) needs to 

be fulfilled only when V 0 is positive, the parametric equation model of Eq. 1 can be further 

refined as: 

20 V eff = V Q - i|(F 0 - V - 6 - A) 4- ^o-^-^) 2 +4^ 0 +2V^A + 2V^ r *A + A 2 |, (12) 

which is the desired enchanced continuous parametric model. 

Using the resulting continuous parametric model of Eq. (12) along with the compensation 
function 9, defined in Eq. (7), a desirable output is achieved whether the input (in the example, 
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V 0 ) is positive or negative. Even in the case when V 0 = 0, the enhanced parametric model of Eq. 
(12) produces the desirable smooth result, which is illustrated in the graph of Fig. 7. 

While the continuous parametric model of Eq. (12) represents a voltage, the terms V efr , 
V 0 , and V can be substituted with other parametric variables such as current (I), Power (P), 
Resistance (R), capacitance (C) and the like. In addition, while the methodology of the present 
invention was demonstrated in connection with the base model of Eq. (1), numerous other base 
models can be transformed into continuous parametric models in accordance with the present 
invention. 

Although the present invention has been described with several embodiments, various 
changes and modifications may be suggested to one skilled in the art. It is intended that the 
present invention encompass such changes and modifications as fall within the scope of the 
appended claims. 
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